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ABSTRACT 

> 

q\ '. We consider the constrained minimal supersymmetric standard model which emerges from 
one theory with a small deviation from Yukawa unification which is adequate for fi > 0. 
We show that this model possesses a wide and natural range of parameters which is 
consistent with the data on b — > sj, the muon anomalous magnetic moment, the cold 
dark matter abundance in the universe, and the Higgs boson masses. 

a: 

Qh' 

o> 1 Introduction 
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We study the phenomenological consequences of imposing on the constrained MSSM 
(CMSSM) an asymptotic relation for the Yukawa couplings at the GUT scale. This 
assumption (Yukawa unification) naturally restricts [1] the t-quark mass to large values 
compatible with the data. Also, the emerging model is highly predictive [2]. Despite of 
its appealing, the simple scheme of a single Yukawa for the three third generations at the 
GUT scale leads to an unacceptable 6-quark mass. This fact excludes minimal versions 
of GUT groups with this property, such as Pati-Salam unification (Gps = SU(4) C x 
SU(2) L x SU(2) R ), 50(10) or E e . 

We consider the SUSY GUT model described in Ref. [3] which is based in the Gps as 
described in Refs.[4] and establish an 'asymptotic' relation for the Yukawa couplings that 
depends on a single complex parameter c: 

h t : h b : h T = |1 + c| : |1 - c\ : |1 + 3c|, (1) 

For simplicity, we will restrict our analysis to real values of < c < 1. The relative 
splitting of the Yukawa couplings becomes: 8h = —(hb—ht)/h t = {h r — h t )/h t = 2c/(l+c). 
This means that the bottom and tau Yukawa couplings split from the top Yukawa coupling 
by the same amount but in opposite directions, with hf, becoming smaller than ht- 
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2 The MSSM with Quasi- Yukawa Unification 



This model, below Mqut, reduces to the MSSM supplemented by the 'asymptotic' Yukawa 
coupling quasi-unification condition in Eq.(l). We will assume universal soft SUSY break- 
ing terms at Mqut, i-e., a common mass for all scalar fields m , a common gaugino mass 
M1/2 and a common trilinear scalar coupling A . In the present work, we will concentrate 
on the fx > 0. The case fx < is phenomenologically less interesting, it will be presented in 
[5] . We follow the notation as well as the RG and radiative electroweak breaking analysis 
of Ref.[6] for the CMSSM with the improvements of Refs.[7, 3] (recall that the sign of fx 
in Refs.[6, 7] is opposite to Ref.[3], which is the one adopted here). 

For any given m b (M z ) in its 95% c.l. range (2.684 - 3.092 GeV for a s (M z ) = 0.1185), 
we can determine the parameters c and tan/3 at M SUSY = (m^m^) 1 / 2 (t 12 are the stop 
mass eigenstates) so that the 'asymptotic' condition in Eq.(l) is satisfied. We use fixed 
values for the running top quark mass m t (m t ) = 166 GeV and the running tau lepton 
mass m T {M z ) = 1.746 GeV and incorporate not only the SUSY correction to the bottom 
quark mass but also the SUSY threshold correction to m T (M S usY) from the approximate 
formula of Ref . [8] . After imposing the conditions of gauge coupling unification, successful 
electroweak breaking and Yukawa quasi-unification in Eq.(l), we are left with three free 
input parameters mo, Mi/ 2 and Aq. In order to make the notation physically more 
transparent, we replace m and Mi/ 2 equivalently by the mass rriLSP (or m^) of the 
lightest supersymmetric particle (LSP), which turns out to be the lightest neutralino (x), 
and the relative mass splitting A f2 = (m f2 — m^)/m^ between the lightest stau mass 
eigenstate (f 2 ) and the LSP. In Fig.l we display the changes on M SUSY and the mass 
of the pseudo-scalar Higgs, m^, for several values of A f2 , and m^iMz)- These changes 
will help us to understand the corresponding predictions for Q LSP h 2 in the presence of 
resonant annihilation channels for values of m A 2 • m LS p. 




Figure 1: The mass parameters tua and Msusy as functions of itilsp, with a s (Mz) — 0.1185 and 
Ao = 0. On the panel of the left we display various values of Af 2 (indicated on the curves) for rrib(Mz) = 
2.888 GeV. On the right panel we fix Af 2 — 1 and show the curves for m\,[Mz) = 2.684 GeV (dashed 
lines), 3.092 GeV (dotted lines) or 2.888 GeV (solid lines). 
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Figure 2: The LSP relic abundance Qlsp h 2 versus m^sp for various Af 2 's (indicated on the curves) 
and with A = 0, m b (M z ) = 2.888 GeV, a s (M z ) = 0.1185. The solid lines (crosses) are obtained by 
micrOMEGAs (our alternative method). The upper bound on £Ilsp h 2 (=0.22) is also depicted. 

Here, the LSP (%) is an almost pure bino. Its relic abundance will be calculated by 
micrOMEGAs [9], which is the most complete code available. It includes all the coannihi- 
lations [10] of neutralinos, charginos, sleptons, squarks and gluinos since it incorporates 
automatically all possible channels by using COMPHEP [11] (A similar calculation has ap- 
peared in Ref.[12].) Also, poles and thresholds are properly handled and one-loop QCD 
corrected Higgs decay widths [13] are used, which is the main improvement provided by 
Ref.[9]. The SUSY corrections [14] to these widths are, however, not included. Fortu- 
nately, in our case, their effect is much smaller than that of the QCD corrections. From 
the recent results of DASI [15], one finds that the 95% c.l. range of VL C dm h 2 is 0.06 — 0.22. 
Therefore, we require that Vl LS p h 2 does not exceed 0.22. 

In order to have an independent check of micrOMEGAs, we also use the following alternative 
method for calculating VLlsp h 2 in our model. In most of the parameter space where 
coannihilations are unimportant, Qlsp h 2 can be calculated by using DarkSUSY [16] 1 . 
Its neutralino annihilation part is in excellent numerical agreement with the recent exact 
analytic calculation of Ref. [18] , the main defect of its current version is that it uses the tree- 
level Higgs decay widths. This can be approximately corrected if, in evaluating the Higgs 
decay widths, we replace m^m^) by m b at the mass of the appropriate Higgs boson in the 
couplings of the 6-quark to the Higgs bosons (see Ref. [9]). In the region of the parameter 
space where coannihilations come into play, the next-to-lightest supersymmetric particle 
(NLSP) turns out to be the T2 and the only relevant coannihilations are the bino-stau 
ones [6, 19] (we do not find bino-stop cannihilations [20], potencially important when A is 
negative). In this region, which is given by A f2 < 0.25, we calculate Qlsp h 2 by using an 
improved version of the analysis of Ref. [6, 7, 21, 23]. The list of bino-stau coannihilation 
channels appropriate for all tan /3's given Ref. [6] has been completed with some additional 
channels as described in [3](see also Refs.[19, 22]). Their corresponding cross sections are 

1 An updated version of this code is now available [17]. 
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Figure 3: Restrictions on the ttilsp ~ Af 2 plane for A = 0, a s (M z ) = 0.1185. From left to right, the 
dashed (dotted) lines depict the lower bounds on tolsp from rriA > 110 GeV, BR(6 — > sj) > 1.9 x 10~ 4 
and Sa^ < 58 x 10~ 10 , and the upper bound on m LS p from Qlsp h 2 < 0.22 for m b (M z ) = 2.684 GeV 
(3.092 GeV). The left (right) solid line depicts the lower (upper) bound on ttilsp from mt > 114.1 GeV 
(Qlsp h 2 < 0.22) for m b (M z ) = 2.888 GeV. The allowed area for m b (M z ) = 2.888 GeV is shaded. 

combined with the results of DarkSUSY as described in [3]. The results presented in Fig. 2 
show an impressive agreement of the two methods. 

We calculate BR(6 — > S7) using the formalism of Ref.[24], where the SM contribution is 
factorized out. This contribution includes the next-to-leading order (NLO) QCD and the 
leading order (LO) QED corrections. The charged Higgs boson contribution to BR(5 — > 
57) is evaluated by including the NLO QCD corrections from Ref.[25]. The dominant 
SUSY contribution includes the NLO QCD corrections from Ref. [26], which hold for large 
tan p. With the considerations of [3] the 95% c.l. range of this branching ratio then turns 
out to be about (1.9 - 4.6) x 10~ 4 . 

According with the latest measurement [27] of the anomalous magnetic moment of the 
muon = (g^ — 2)/2, the deviation of from its predicted value in the SM [28], 5a^, 
is found to lie, at 95% c.l., in the range from -4.7 x 1(T 10 to 56 x 1(T 10 when the SM 
calculations based in e + e~ data and in r data are both taken into account. The calculation 
of 5 dp in the CMSSM is performed here by using the analysis of Ref. [29], the updating 
of the experimental bounds does not introduce significant differences respect the results 
presented in [3]. 

We will also impose the 95% c.l. LEP bound on the lightest CP-even neutral Higgs boson 
mass m h > 114.1 GeV. In the CMSSM, this bound holds almost always for all tank's, 
at least as long as CP is conserved. The CP-even neutral Higgs boson mass matrix by 
using FeynHiggsFast [30]. Finally, for the values of tan/3 which appear here (about 60), 
the CDF results yield the 95% c.l. bound m A > 110 GeV [31]. 
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3 The Allowed Parameter Space 

The restrictions on the rriLSP — A ?2 plane, for A = and with the central value of 
a s (M z ) = 0.1185, are shown in Fig.3. The dashed (dotted) lines correspond to the 95% 
c.l. lower (upper) experimental bound on m b (M z ) which is 2.684 GeV (3.092 GeV), while 
the solid lines correspond to the central experimental value of mb(M z ) = 2.888 GeV. 
From left to right, the dashed (dotted) lines depict the lower bounds on itllsp from the 
constraints m A > 110 GeV, BR(b -> sj) > 1.9 x 10~ 4 and 5a^ < 58 x 10~ 10 , and the 95% 
c.l. upper bound on m L sp from Qlsp h 2 < 0.22. The constraints BR(6 — > 57) < 4.6 x 10~ 4 
and 5a^ > — 6 x 10~ 10 do not restrict the parameters since they are always satisfied for 
pi > 0. The left solid line depicts the lower bound on m LS p from m h > 114.1 GeV which 
does not depend much on m b (M z ), while the right solid line corresponds to Qlsp h 2 = 0.22 
for the central value of m\ ) {M z ). We see that is always smaller than Ivtilsp but close 
to it. Thus, generally, the neutralino annihilation via the s-channel exchange of an A- 
boson is by far the dominant (co) annihilation process. We also observe that, as itilsp 
or Af 2 increase, we move away from the A-pole, which thus becomes less efficient. As a 
consequence, Qlsp h 2 increases with uilsp or Af 2 (see Fig. 2). 

In the allowed (shaded) area of Fig. 3 which corresponds to the central value of rrib(M z ), 
the parameter c (tan/5) varies between about 0.15 and 0.20 (58 and 59). For rrib(M z ) 
fixed to its lower or upper bound, we find that, in the corresponding allowed area, the 
parameter c (tan/5) ranges between about 0.17 and 0.23 (59 and 61) or 0.13 and 0.17 
(56 and 58). We observe that, as we increase rrib(M z ), the parameter c decreases and 
we get closer to exact Yukawa unification. This behavior is certainly consistent with the 
fact that the value of rrib(M z ) which corresponds to exact Yukawa unification lies well 
above its 95% c.l. range. The LSP mass is restricted to be higher than about 138 GeV 
for A = and a s (M z ) = 0.1185, with the minimum being practically A f2 -independent. 
At this minimum, c « 0.16 - 0.20 (c « 0.13 - 0.23) and tan/5 « 59 (tan/5 w 58 - 61) for 
m b (M z ) = 2.888 GeV (m b (M z ) = 2.684 - 3.092 GeV). 

In Fig. 3, we present the restrictions on the ttilsp — A /Mi/ 2 plane for m b (M z ) = 
2.888 GeV, a s (M z ) = 0.1185 and fixed values of A f2 . The most significant restriction on 
the allowed area are due to the displacement of the Qlsp h 2 = 0.22 line as A ?2 increases, 
showing clearly the effect of the bino-stau coannihilations on the left panel. On the right 
panel we can observe that the allowed area becomes narrower as |A /Mi/ 2 | 7^ 0. 

4 Conclusions 

We showed that, in the particular model with Yukawa quasi- unification considered, there 
exists a wide and natural range of CMSSM parameters which is consistent with all the 
above constraints. We found that, within the investigated part of the overall allowed pa- 
rameter space, the parameter tan/5 ranges between about 58 and 61 and the 'asymptotic' 
splitting between the bottom (or tau) and the top Yukawa couplings varies in the range 
26 — 35% for central values of rrib(M z ) and a s (M z ). 
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Figure 4: Restrictions on the m L sp - A /M 1/2 plane for Af 2 = (left) and 1 (right), m b (M z ) = 
2.888 GeV, a s (Mz) — 0.1185. The solid, dashed, dot-dashed and dotted lines correspond to the lower 
bounds on m LSP from m A > 110 GeV, BR(b -> s 7 ) > 1.9 x 10~ 4 , 5a^ < 58 x 10~ 10 and m h > 114.1 GeV 
respectively. The upper bound on itilsp from Qlsp h 2 < 0.22 does not appear in the left panel since 
it lies at itilsp > 500 GeV. The allowed area is shaded. For comparison, we also display on the left 
panel the bounds from Qlsp h 2 < 0.22 (double dot-dashed lines) for A ?2 = 0.1 and 0.03, as indicated. 
The upper bound on m LSP from the cosmological constraint Qlsp h 2 < 0.22 corresponds to the double 
dot-dashed line on the left panel. 
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